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Abstract 

In this paper we derive an integral equation for the evolution of unin- 
tegrated (longitudinally) polarized quark and gluon parton distributions. 
The conventional CCFM framework is modified at small x in order to 
incorporate the QCD expectations concerning the double ln 2 (l/x) resum- 
mation at low x for the integrated distributions. Complete Altarelli-Parisi 
splitting functions are included, that makes the formalism compatible with 
the LO Altarelli-Parisi evolution at large and moderately small values of x. 
The obtained modified polarized CCFM equation is shown to be partially 
diagonalized by the Fourier-Bessel transformation. Results of the numeri- 
cal solution for this modifed polarized CCFM equation for the non-singlet 
quark distributions are presented. 



1 Introduction 



The basic, universal quantities which describe the inclusive cross-sections of hard 
processes within the QCD improved parton model are the scale dependent par- 
ton distributions. These distributions depend upon the longitudinal momentum 
fraction x and the hard scale Q 2 and correspond to the integrals over transverse 
momentum k± of the so called unintegrated distributions describing the (scale 
dependent) x and k± parton distributions. The unintegrated distributions are 
needed in the description of less inclusive quantities which are sensitive to the 
transverse momenta of the partons [1] - [5]. 

The unintegrated parton distributions are described in perturbative QCD by 
the Catani, Ciafaloni, Fiorani, Marchesini (CCFM) equation [6, 7, 8] which is 
based on color coherence which implies angular ordering. It embodies in a uni- 
fied way the (LO) Altarelli-Parisi evolution at large and moderately small values 
of x with the BFKL dynamics at small x. 

The CCFM equation which was originally formulated for unpolarized parton 
densities has recently been generalized to spin dependent unintegrated parton 
distributions [9]. Novel feature of the spin dependent case is its different structure 
at small x and in particular the absence of the non-eikonal (or non-Sudakov) form- 
factors. We include here also the (spin dependent) quark distributions, leading 
first to a system of CCFM equations for unintegrated spin dependent quark and 
gluon distributions. We shall then extend the analysis of Ref. [9] along the 
following lines: 

1. We shall incorporate theoretical QCD expectations concerning the double 
logarithmic small x effects. We observe in particular that both the angu- 
lar ordering constraint and the kernels of the corresponding system of the 
CCFM equations have to be appropriately modified. 

2. We shall include complete splitting functions P a b{z) and not only their 
singular and finite parts in the limit z — > 1 and z — > respectively. 

3. We shall explore the relative simplicity of the system of the modified polar- 
ized CCFM equations due to absence of non-eikonal form-factors and utilize 
the transverse coordinate representation of those equations. The transverse 
coordinate b± is related to the transverse momentum of the partons through 
the Fourier-Bessel transformation and the system of the CCFM equation 
can be partially diagonalized by this transformation. 

The content of our paper is as follows: in the next section we recall the original 
formulation of the CCFM equation developed in [9] and include the polarized 
quark distributions, while in sections 3 and 4 we formulate the modifications of the 
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CCFM framework, which will incorporate the QCD expectations concerning the 
double ln 2 (l/x) resummation at low x and the complete Altarelli-Parisi evolution 
of integrated densities in moderately small and large values of x. We observe that 
the modifed CCFM equations, including their extension discussed in sections 4 
and 5 can be partially diagonalized by the Fourier-Bessel transformation. Section 
5 is devoted to the numerical analysis of the modified CCFM equations and for 
simplicity we limit ourselves to the non-singlet quark distributions. Finally, in 
Section 6 we summarize our main results and give our conclusions. 



2 The original CCFM formulation for longitu- 
dinally polarized unintegrated parton distri- 
butions 

Along the lines of [6] [7] [8] a version of the CCFM equation has been derived 
for the longitudinally polarized gluon distribution [9]. For the corresponding 
expressions including quarks one has to take into consideration that the soft 
emission occurs predominantly from a ladder of gluons because of their larger 
spin quantum number. Therefore, quarks appear only as initial states playing 
automatically the role of the hardest emission which enters in the Altarelli Parisi 
splitting function, see Fig. 1. Therefore, the arguments presented in [9] are valid 
in the case of quarks as well. The polarization does not enter into the soft 
emission. Consequently, again a non-eikonal form factor is absent and the eikonal 
form factor is identical to the one in the unpolarized case. Furthermore, in its 
original formulation the CCFM equation is an evolution equation valid only in 
the limits x — > and x — > 1, and therefore, the Altarelli Parisi kernels enter 
only in interpolated form, taking the limits z — > and z — > 1 into account. The 
complete polarized CCFM equations take the following general form: 



Af(x,kl,Q 2 ) = (AP CCFM ®Af)(:r,A;i,Q 2 ) 
Af k (x,k\,Q 2 ) = (ApC q CFM ®Af k )(x,klQ 2 ) . (1) 

Here we have used: 



/ A/ x /AP™ AP™\ 



Af 



and AP CCFM 

U/J Vap™ AP™y 



(2) 



Here Af g is the unintegrated polarized gluon distribution function, A/s the un- 
integrated polarized singlet quark parton distribution function and Afk denotes 
any polarized non-singlet combination of the quark parton distribution functions 
like the triplet (A/ 3 ) and the octet contribution (A/ 8 ). For the convolution one 
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has the structure: 



(AP CCFM 8 Af )( x , kl, Q 2 ) = f" d -^r/4 -6(0 - 4<ul 

JO Z7T JQi Q I Jx z 



r2n d9 q± r°° dqj_ f 1 -Q»/idz i 
2ir JqI q\ 
xAP CCFM (z,gi,Q 2 )A/(^,^ 2 ,gi), (3) 



using k' ± — k± + (1 — z)q±. In terms of the principles discussed above one obtains 
for the CCFM splitting kernels: 

AP™(z,q 2 ± ,Q 2 ) = A(*>(Q 2 , {z q lf f Ml{ l- Z)2) AP^(Z) 
AP™™(z, q 2 ± ,Q 2 ) = A^>(Q 2 , {zqlf) as{ql{l 2 - Z)2) AP^\z) 
AP™(.,gi,Q 2 ) = A^(Q 2 , (zqlf f'^- Z?) AP^{z) 
AP™(,,, 2 ,Q 2 ) = A^Q 2 , {zqlf) aMl{l 2 - Z)2) AP q f\z) . (4) 
The eikonal form factors have the form: 

,/2 



{zq) ) = exp L IF I — * ) {S) 

and the interpolating Altarelli-Parisi kernels read then (c. f. [10]): 

AP^\z) = 



ap 9 a /°w = i* 



1-2 
r 



1 + 2 2 

-^j = ( /— — 

jAPO 



AR A /°(z) = C 



2 

1 - 

F 



AP g A r(z) = C F (2-z). (6) 

The fact that in the polarized case the CCFM splitting kernels are independent 
of k± allows a factorization in terms of the Fourier-Bessel transformation: 

(AP CCFM ® Af)(x,b 2 ± ,Q 2 ) = jd 2 k ± e~* b ^(AP CCFM ®Af)(x,kl,Q 2 ) 

fidz rQ 2 /^ dq 2 ±T ... ... 1 1 h 

Jx z jq^ qj_ 

xAP CCFM (z,qlQ 2 )Af(x/z,bl,ql) . (7) 

One should note that f(x,Q 2 ) = f(x,0,Q 2 ) is just the k± integrated parton 
distribution. 
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Figure 1: Construction of the polarized CCFM equation for quarks from the one 
for gluons. The kernels for qq, qg and gq splitting are subsequently constructed 
from the gg soft emission by amending the corresponding hard initial states. Here 
the solid line describes gluons while the dashed line denotes quarks. 
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3 Modifications to make contact with ladder di- 



agrams 

It may easily be observed that the angular ordering constraint which is embodied 
within the CCFM equation generates the double ln 2 (l/:r) terms for the inte- 
grated parton distributions. The result concerning resummation of those terms 
turns out to be, however, different from the QCD expectations discussed in [11], 
[12], [13], [14]. The double ln 2 {l/x) effects in QCD are generated by the ladder 
and non-ladder bremsstrahlung diagrams. 

In order to get the expected double logarithmic limit of the integrated distri- 
butions corresponding to the ladder diagrams contribution it is sufficient to re- 
place the angular ordering constraint Q(Q — z\q±\) by the stronger constraint 
0(Q 2 — zq\) in the corresponding evolution equations for integrated distribu- 
tions. The latter just correspond to the Fourier transformed unintegrated evo- 
lution equations at bj_ = 0. We assume that this replacement can be done for 
arbitrary values of the transverse coordinate b±. Due to this replacement we will 
obtain an equation which incorporates the known collinear (LO Altarelli-Parisi) 
evolution for a not too small x and the double logarithmic asymptotics for x <C 1. 
However, the equation does not sum up the single ln(l/a;) contributions as it was 
done by the original CCFM/BFKL (unpolarized) equation. Having this point in 
mind we will henceforward call the evolution equation for the polarized uninte- 
grated parton distributions which we derive starting from the CCFM formulation 
and where we are including our modfications a modified polarized CCFM evolu- 
tion equation. Apart from this substitution in the argument of the theta function 
to include the double logarithmic contributions, we shall also make the following 
modifications of the original CCFM equations proposed in [9]: 

1. The argument of a s will be set equal to q\ instead of q\(l — z) 2 . 

2. The non-singular parts of the splitting function(s) will be included in the 
definition of the Sudakov form-factor (s). 

3. Following Ref. [14] we include complete splitting functions P a b(z) and not 
only their singular parts at z — 1 and constant contributions at z = 0. 

4. We represent the splitting functions AP ab (z) as AP ab (z) = AP ab (0) + 
AP ab {z) where AP ab (0) = 0. Following [14] we shall multiply AP ab (0) 
and AP ab (z) by Q(Q 2 — zq\) and Q(Q 2 — q±) respectively in the integrands 
of the corresponding integral equations. Following the terminology of Ref. 
[14] we call the corresponding contributions to the evolution kernels the 
'ladder' and 'Altarelli - Parisi ' contributions respectively. 

5. We shall 'unfold' the eikonal form factors in order to treat real emission 
and virtual correction terms on equal footing. 
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Using those prescriptions we get the following unfolded and modified polarized 
CCFM equations: 



Af g (x,k 2 ± ,Q 2 ) = Af°(x,k] 



d 2 q_L a s 



-e(9i-0g) / 

J 



1 dz 

z 



x< Q(z - x) 



(l20(Q 2 - zql) + 60(Q 2 - ql) - 2zjj Af g (x/z, k'l, q\) 



+ [1&(Q 2 - zq 2 ± ) - -z&(Q 2 - ql)) AMx/z, k'l, ql) 
-zQ(Q 2 - ql) (JL - y + f ) A/ fl (x, kl, ql)\ 



(8) 



A/ E (x, fci, Q 2 ) = A/°(x, fci) + / ^t^O(ql - Ql) 



d 2 q± a s (ql), 

Tiq 2 27T 



1 dz 



z 



x< 6(2 - x) 



-N F Q(Q - ^i) + 2zN F Q(Q 2 - q 2 ± ))Af g (x/z, k'l, ql) 



+1 (e(Q 2 - *£) + ^e(g 2 - <£) ) a/ s (x/z, ^ 2 , ql) 



-zQ{Q 2 -ql) 



3 1-2 



2 A/ E (x,fci,gl) , 



(9) 



where we have put explicit numbers for the factors Ca and Cp, and also intro- 
duced the singlet spin dependent unintegrated quark distributions: 



A/ E = E&(A/ B + A/ & ) 



(10) 



The 'non-singlet' quark distributions evolve as A/ s but without the gluon con- 
tribution on the r.h.s. of the integral equation. 

In equations (8) and (9) we set the upper limit of integration over dz equal 
to 1 instead of 1 — Qo/q since the integrands are free from singularities at z — 1. 

— * — * 

It should be noted that k' ± = k± + (1 — z)q±. It should also be noted that the 
inhomogeneous terms Af g and A/£ in equations (8) and (9) do not contain the 
Sudakov form-factors. They can be chosen to have (for instance) the Gaussian 
form in k± normalized to unity and multiplied by the input (integrated) distri- 
butions at the scale Q 2 = Qfo. The latter could be taken from one of the existing 
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(LO) QCD analysis of spin dependent parton distributions. To be precise the in- 
homogeneous terms are related to the starting distributions at the reference scale 
Ql in the 'single loop' approximation [15, 16] corresponding to the replacement 
Q(Q 2 — zql) by Q(Q 2 — ql) in the integrals in equations (8,9), since in general 
the integrals containing the function <d(Q 2 — zql) do not vanish at Q 2 = Q%. 
Parameterization of the driving term in terms of the parton distribution may be 
regarded as a reasonable approximation, particularly in the region of large and 
moderately small values of x which is dominated by the single loop approximation. 

Taking the Fourier-Bessel transformation on both sides of equations (8) and 
(9) we get the following equations for the distributions f g (x, b 2 ± , Q 2 ) and /s(x, b 2 ± , Q 2 ): 



f 9 (x, b 2 ± , Q 2 ) = f g {x, b 2 ± ) + f ^SMLq^ ± - Ql) f f 

J 7TQ I 2,71 JO Z 



- 1 dz 



irql 2tx 



x\j [b ± (l-z)q ± }e(z-x) 



(l2Q(Q 2 - zql) + 60(Q 2 - q\) - 2zj) f g (x/z, b 



',2 



+ (J0(Q 2 - zq 2 ± ) - ^zQ(Q 2 - ql) j fc(x/z, b 2 ± , ql) 
-zQ(Q 2 - ?i) - y + f ) /„(*, bl, q 2 ± ) \ 



(11) 



h(x, b 2 ± , Q 2 ) = M(x, b 2 ± ) + / ^^Q(ql - Ql) I' - 

J q i Air Jo z 



- 1 dz 



ql 2tt 



x{j [b ± (l-z)q ± ]e(z-x) 



N F Q(Q 2 -zql)+2zN F e(Q 2 -ql))f g (x/z, b 2 ± , ql) 



+ 3 [0(Q 2 - zql) + YTJ®(Q 2 - Ql) Mx/z, bl, ql) 



-zB(Q 2 - q 2 ± ) 



3(1 -z) 



2 )Mx,bl,ql) , 



where: 



(12) 



Ux,bl,Q 2 ) = j d 2 k ± exp(-tk ± -b ± ) Af t (x,kl,Q 2 ) 

POO 

= 2tt / k ± dk ± J (k ± b ± )Af t (x,kl,Q 2 ) , 
Jo 



(13) 



and 



A/ifofcLQ 2 ) = J ^y 2 exp(tk ± -b ± )f t (x,bl,Q 2 ) 



(27T) 



1 z 100 
2n Jo 



r 'b ± db ± J {k ± b ± )Afi{x,bl,Q 2 ) . (14) 

JO 



From a physical point of view one can interpret b± as an 'impact parameter' giving 
the transverse distance of the partonic probe. Then the following expressions are 
obtained: 

J g (x,b 2 ± ,Q 2 )=T g (x,b 2 ± ) 

r 1 dz [Q 2 / z dq 2 , a s (q 2 ,)- . . , 9 ,.8 T ... , xx 

+ / , T/r"^ S # i 7 f (^,H,i)12Jb(|bJ|q,|(l-,)) 
ix z JQl q ± An 

(ladder) 



*dz fQ 2 dq 2 ± a s (ql). 



+ 



JiL 

J x Z J Q 

r 1 (h rQ 2 dq 2 ± a s (q 2 ± ) 

J x Z J C 



+ L ~l Ql ^r^f i 7=(*A.ft.d)(-^)*(|bx||q 1 |(i-*)) 



Ql Q± 2tt 



Qz 



f g (x/z,b 2 ± ,q 2 ± ) -f g (x,b 2 ± ,q 2 ± ) 
(1-z) 

-2j(x/z,b 2 ± ,ql) 



J (|bJ|qJ(l-^)) 



Q 2 dql a s (ql) 
Ql q± 2tt 



f g (x,b 2 ± ,q 2 ± ) 



11 N f 



+ 61n(l -x) 



(Altarelli Parisi) 

J^x ) b\ ) Q 2 )=%{x ) bl) 
^dz [Q 2 /* dq 2 ± a s {q 2 ± )- 



(15) 



z JqI q 2 ± 2tt 



^dz rQ 2 /-dq\ a s (q] 



x z JqI q]_ 2n 



f^x/z,bi,qi)-J (\b ± \\ q± \(l-z)) 



f g (x/z,bi,ql)N F J (\b ± \\q ± \(l-z)) 
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(ladder) 



+ 



dz fQ 2 dq\ a s (q 2 A 



z JqI qi 2tt 



■Jo(|b ± ||q ± |(l-z)) 



x • 



(z + z 2 )f^(x/z, b\, q\) - 2zfx(x, b 2 ± , q\) 



f 



:i - z) 



2 + -\n(l-x) 
o 



+ — 



1 dz /•(-'" dq' J <\Jq 2 ) - f , ;2 , 
z JqI q 2 ± 2n lg 



f Q (x/z,bl, q l)2zN F J (\b ± \\q ± \(l-z)) 



(Altarelli Parisi) 



(16) 



The contribution for the quark non-singlet part can be simply obtained from the 
expressions for the singlet part leaving simply out all gluonic contributions: 

fq Ns( X i b 2 ±, Q 2 ) = f q Ns( X ' &J_) 

*dz fQ 2 /* dqla s {q 2 ^ A 



+ 



fil 

J x Z JL 



faNs(x/z,b ± ,q ± )- J (|b ± ||qJ(l - z)) 



Ql q± 2tt .» 



(ladder) 



+ 



idz fQ 2 dql a s (ql) 
Ql Q 2 ± 2tt 



fit 

J x Z J L 



Jo(|bx||cLi|(l-^)) 



4 
X 3 



{z + z 2 )f q NS {x/z, b 2 ± , ql) - 2zf q NS (x, b\, ql) 



+ 



Q 2 dql a s (ql) 
f Q± 2 ^ 



fq Ns( X i &J_> Q±) 



2 + -ln(l -x) 

o 



(Altarelli Parisi) 



(17) 



The expressions show that except for the scale of a s and the occurrence of 
the Bessel function Jo the expression match exactly the ones derived in [14] 
calculating ladder diagrams and combining this with the Altarelli Parisi evolution. 
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This shows the tight relationship that can be shown between modified CCFM and 
the ladder contributions. 



4 Inclusion of non-ladder diagrams 

There is a third contribution to the evolution of unintegrated parton distribu- 
tions which is not covered by the 'Altarelli - Parisi +ladder' approximation of 
the modified polarized CCFM equation, these are the non-ladder bremsstrahlung 
contributions. The method of implementing the non-ladder bremsstrahlung cor- 
rections in general into the double logarithmic resummation was proposed by 
Kirschner and Lipatov [13]. For integrated polarized parton distributions they 
have been implemented in Ref. [14]. The method developed in Ref. [13] is based 
on the infrared equations for the partial waves F , 8 (cj, a s ): 



Fo,s 



The infrared equations for the partial waves read 



(18) 



Ana* 

47TQU 



-Mr 



2a, 



-M 8 + 



TIUJ 

a.N d 



2 F 8 (w,a s )G + 
F 8 (^,a s ) + 



8tt 2 cj 
1 



Fo(^,a s ) 



u> 2ituj duj 

where the matrices M , M 8 and Go are given by: 



Fo(cu, a s ) 



(19) 
(20) 



-l 

2N 



N F 
2 



N 2N 



(21) 



/ AP qq (0) AP qg (0) 



M = 



\AP gq (0) AP gg (0) 




G = 



/ N 2 -l 
2N 



\ N 



(22) 



(23) 



The singlet partial wave matrix Fo is linked with the anomalous dimension ma- 
trix / -fg ES (uj, a s ) controlling the evolution of the moments of the integrated spin 
dependent distributions: 

2„.RES ( , . „, ^ (24) 



Fo = 87T i 7 ^(^a s ) 
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The anomalous dimension matrix corresponding to the solution of equation (19) 
is given by: 



7 f B5 (a;,a.) = =-ll- 



\ 



2a s /M F 8 (u,a s )G ^ 



1TLU 



2ll 2 LU 2 



(25) 



The anomalous dimension matrix contains the resummation of the double loga- 
rithmic ln 2 (l/x) effects which correspond to the sum of powers of a s /co 2 in the u 
space. 

It should be noted that the inhomogeneous term in the nonlinear equation (19) 
which is proportional to Mo is also proportional to the kernel matrix defining 
the ladder diagram contributions in the double logarithmic approximation to the 
evolution equation for the (integrated) parton distributions (cf. equations (8, 9) 
at b± = 0). The fact that in equation (19) M appears in the inhomogeneous 
term while in equations (8, 9) it appears in the kernel matrix is linked with the 
fact that equation (19) defines the gg-scattering amplitude in the u representation 
while equations (8, 9) define the parton distributions. The double logarithmic 
approximation of equations (8, 9) at b±_ = would generate the anomalous di- 
mension matrix Tedder corresponding to ladder diagrams in the double logarithmic 
approximation which is given by: 



l^L(uj,a a ) = - (1 



2a s M 



tyuj z 



(26) 



It has been observed in Ref. [14] that in order to get complete anomalous dimen- 
sion given by equation (25) one has to add the corresponding terms proportional 
to F 8 (a;,a s )Go in the kernel matrix defining contribution of ladder diagrams. 
In the two-scale unintegrated case one can simply add them by analogy to the 
Altarelli - Parisi and ladder contribution by inserting the factor Jo(b±q±(l — z)). 
The results are given below: 

7Mbl,Q 2 )=fjx,bl) 



1 dz rQ 2 / z do 2 , a s (q 2 )- , , 9 9n 8 t/1 , 
; _j _^_^±l f ^ b 2_^ J (| b± ||qJ(l-z)) 



+ 



z JqI q\ 2tt 
idz [Q 2 /* dqlasiql 



fil 

Jx Z JL 



Ql q\ 2n '" 
(ladder) 



fJx/z,bi,qi)12J (\b ± \\ q± \(l - z)) 



+ 



l dz [Q 2 dq 2 ± a s (q 2 ± )- 



2tt 



f s (x/z,bl,q 2 ± ) - -z J (|b ± ||q ± |(l-z)) 
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+ f'dz rQ 2 dqj as(ql) 6z 
Jx z JqI qj_ 2n 



f q (x/z,bl,q 2 ± ) - f Jx,b 2 ± ,q 2 ± ) 



(1-z) 



-2fJx/z,bi,qi) 



Jo(|b ± ||<u|(l-^)) 



[Q 2 dq\ a s (qj) - 2 2 
(Altarelli Parisi) 



11 iV P 



+ 61n(l -x) 



" l dz fQ 2 dqla s (q 



-JtL 

Jx Z JO, 

r 1 dz r 

Jx Z JQ° 



Ql q 2 ± 2tt 
l dz fQ 2 / z dq\ a s (qi) 



J (|bJ|qJ(l-z))^ 
J (|b ± ||q ± |(l-z)) 



■-)^r) f^x/z^ql) 



2tt 2 



91 



U 2 



x z Jq 2 qj_ 2-k 
(non — ladder) 



fz(x,bl,Q 2 ) = fl(x,b 2 ± 



- 1 dz fQ 2 / z dq 2 a s (q 2 A- . . , 9 o x 4 T .„ . 
- / — / . -^^^/ s (x/z,6i, g i)-J (|b ± ||q ± |(l - z)) 



Q 2 ) 2tt 2 I Jg 



f(x/z,b 2 ± ,ql) 



99 



x Z 



qi 2tt 



1 dz [ Q2 / z dq 2 ± a s {q 2 ± )^ 



z Jq 2 qi 2tt 
(ladder) 



/ g (z/z,&i, ? i)AW (|bJ|q ± |(l-z)) 



(27) 



+ r^[ Q2d 4^j ( Mq± \(i-z)) 



z Jq 2 qj_ 2-k 



x- 



z + z 2 )fdx/z, 61, q 2 ± ) - 2zf s (x, 61, q 2 ± ) 



f 

Jc 



(1-z) 



2 + -ln(l -x) 



4 I 1 - / Q2 ^^4^L(^,6l,gi)2ziV F J (|bJ| q J(l-^) 



lx z Jq 2 q 2 ± 2tt 
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(Altarelli Parisi) 



1 dz [Q dq\ a s (q 



-JtL 

J x Z J L 



Ql q\ 2tt 
1 dz [Q 2 / Z dqi a s {q\ 



r L dz r 

J x Z J C 



Q 2 qj_ 2tt 
non — ladder) 



J (|bJ|qJ(l-z))^ 
Jo(|b ± ||q ± |(l-z)) 



u 2 



f^x/z^ql) 



2tt 2 



n 



u 2 



Q 2 J 2n 2 



f g (x/z,b ± ,q ± ) 



19 



(28) 



The contribution for the quark non-singlet part can be simply obtained from the 
expressions for the singlet part leaving simply out all gluonic contributions: 

fq Ns( x i Q 2 ) = f q ns( X ' fr±) 

^dz rQ 2 /*dql a s {q 2 ^ A 



+ 



fit 

J x Z J L 



Ql q\ 2vr ^-^-,-,-^^3 



(ladder) 



f q Ns( x / z , b ±,Q±)o J o(\ h ±\\ ( l±\( 1 ~ z )) 



x z JqI qj_ 2tx 



X- 



z + z )fg ns( x / z , bl, ql) - 2zf q NS (x, b 2 ± , q 2 ± ) 



+ 



? 2 dq 2 ± a s (ql 



'fq Ns( X i ^±5 Q±) 



(1-z) 



2 + -ln(l -x) 
o 



(Altarelli Parisi) 



1 dz fQ 2 dq 2 ± a s (q 



-JtL 

J x Z J L 



Ql ql 2tt 

1 dz [Q 2 /* dq 2 ± a s (ql) 
IT 2tt 



f L dz r 

J x Z J C 



J (|bJ|qJ(l-z))^ 
J (|b ± ||q ± |(l-z)) 



(z) 



Go 



f q Ns( x /z,b ± ,q ± ) 



u 2 



2ir 2 Jq 

/ qq 
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(29) 



Here 




(z) is the inverse Mellin transformation of F 8 (u;)/u; 2 . Derived for 



large iV and fixed a s , one can use an approximate form [11]: 

[z) = 27ia s M 8 \n\z) , 



FBorn 
8 



LO 



(30) 



so that one gets in our case e.g.: 



Co 2 



Go 

2tt 2 



N 2 



qq 



4irN 2 



-a s (qi)\n 2 (z) . 



(31) 



5 Numerical studies 
5.1 The input distributions 

For the polarized input distributions we are going to use the LO Standard Sce- 
nario parameterization given in [17]. The structure of the integrated input dis- 
tributions is a factor multiplied with the unpolarized input distributions given 
in [18]. As we have 'de-exponentiated' the eikonal form factor in the evolution 
equations, the input distributions have to be Q 2 independent. This is different 
from the usual unpolarized CCFM-input distributions, where one uses for the Q 2 
dependence the corresponding eikonal form factor [19]. For the k± dependence of 
the input distributions a Gaussian Ansatz is common [19]. In this way the input 
distribution as defined at an input scale Ql have the following general scheme: 



Af?(x,k 2 ± ,Q 2 ) = A Pl (x,Ql)exp 
Af?(x,b 2 ± ,Q 2 ) = A Pl (x,Ql)exp 



k 2 



a 2 / 7RX 2 



b 2 a 2 " 



(32) 



where A Pi (x, Qq) are the input integrated spin dependent parton distributions at 
the reference scale Q^. 



5.2 Features of the evolution 

Setting b± = in Eqs. (12) and (11), one obtains evolution equations for the inte- 
grated polarized parton distributions equivalent to those given in [14]. Therefore, 
the x and Q 2 dependence of the unintegrated parton distributions will be exactly 
the same as given by the Altarelli Parisi evolution equations supplemented by 
ladder and non-ladder contributions. In this way the modified polarized CCFM 
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equation as discussed here is consistent with the standard evolution of the inte- 
grated polarized parton distributions. In fact the whole transverse momentum 
dependence is governed by the inclusion of the factor Jo(b±q±(l — z)). In principle, 
the only genuine new feature of the modified CCFM equation is the k± (or b±) 
dependence. At large and moderately small values of x one can make the 'single- 
loop' approximation [15, 16] corresponding to the replacement Q(Q 2 — zq\) by 
just Q(Q 2 — g 2 ) in the 'ladder' contribution and to neglecting the 'non-ladder' 
contribution. At = the modified polarized CCFM equation in the single 
loop approximation equation reduces then to the LO Altarelli-Parisi equation in 
the integral form. 

Numerically, we perform the evolution using the Chebyshev approximation 
technique as discussed in the Appendix. The calculation is very time consuming, 
therefore, we take only 8 polynomials in all cases into account. As the Bessel 
function J is oscillating the corresponding integration is numerically quite prob- 
lematic. As a pragmatical solution the integration is only performed up to the 
fourth zero. We have checked that this procedure provides stable results. As 
the mathematical structure is the same for singlet and non-singlet contributions 
we can restrict ourselves for the discussion of the k± dependence to the simple 
non-singlet case. In Fig. 2 we show the evolution of the k± dependence for the 
triplet contribution A/ 3 = | (^Au + Au — Ad — Adj . The input distributions at 

Ql = 0.26 GeV 2 taken from the GRVS LO Standard Scenario set [17] are com- 
pared to the evolved distributions at Q 2 = 10.0 GeV 2 and Q 2 = 100.0 GeV 2 . The 
width of the initial transverse momentum dependence a has been chosen to be 
1 GeV. For the simulation the full content of the equation, the Altarelli-Parisi, 
ladder and non-ladder contributions all have been included. It is seen that due 
to the evolution the k± dependence is broadening away from a Gaussian behavior 
to a more purely exponential decay. Such a feature has already been seen in the 
purely gluonic formulation of the genuine polarized CCFM equation as discussed 
in [9]. 

6 Summary and conclusions 

In this paper we have derived an integral equation of the unintegrated polarized 
parton distributions starting from a genuine polarized CCFM (pCCFM) formu- 
lation and including some modifications incorporating the known collinear (LO 
Altarelli-Parisi) evolution for a not too small x and the double log asymptotics 
for x <C 1. This evolution equation which we call modified polarized CCFM 
equation does not sum up the single ln(l/x) as it was done by the CCFM/BFKL 
unpolarized equation. An inclusion of the single log contributions would be very 
interesting, but is unfortunately beyond the scope of this article. The modified 
CCFM equation yields an approximate description, which contains contributions 
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Figure 2: Transverse momentum dependence for the triplet contribution A/3 
= I (Aw + Au — Ad — Adj of the full modified pCCFM evolution (approximate 
form: Altarelli Parisi + ladder + non-ladder). The thin lines show the input 
distributions GRSV LO Standard Scenario [17] for Q 2 = 0.26 GeV 2 for x = 0.1 
(solid) and x = 0.01 (dashed), while the bold lines show the same distribution 
evolved to Q 2 = 10 GeV 2 (top) and Q 2 = 100 GeV 2 (bottom). It can be seen 
from the logarithmical scale that the evolution leads to a k± broadening away 
from the Gaussian shape to a mere simple exponential decay. 
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Figure 3: Relationship between the pCCFM evolution and its modified approxi- 
mate form for unintegrated parton distributions on the one hand and the evolu- 
tion of integrated parton distributions on the other hand. 
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of the type of 'Altarelli Parisi', 'ladder' and 'non-ladder' and which correspond 
to respective contributions in the integrated case. We have utilized the fact that 
these equations can be diagonalized using the Fourier-Bessel transformation with 
the 'impact parameter' b±. The nice feature of this transformation is that by 
simply setting b± = one obtains already the corresponding expressions for the 
integrated parton distributions. In fact, the difference between both is only a 
factor Jo(b±q±(l — z)). As to the b± dependence the mathematical structure is 
the same for the non-singlet and singlet contributions to the unintegrated par- 
ton distributions. Using the technique of the Chebyshev approximation we have 
performed the evolution for the triplet distribution A/ 3 and observed a char- 
acteristic k± broadening away from the Gaussian form to a more exponential 
decay. The relationship between the genuine and the modified pCCFM for un- 
integrated parton distributions on the one hand and the evolution of integrated 
parton distributions on the other hand is displayed in Fig. 3. The modifications 
to the genuine pCCFM equation correspond to the Altarelli-Parisi part plus lad- 
der contributions in the evolution of the unintegrated parton distributions. The 
non-ladder contributions have to be added by hand. As the only difference in the 
evolution kernels for the integrated and the unintegrated parton distributions in 
Fourier space is the factor Jo(b±q±(l — z)), which becomes unity as the transverse 
impact parameter b± goes to zero, it is clear that the diagram between transverse 
integration and evolution shown in Fig. 3 commutes. 

The next step will be to use the modified polarized CCFM equation as presented 
here to construct a consistent set of unintegrated polarized parton distributions 
from existing data. 
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Appendix 

Approximation by Chebyshev polynomials 

A possible diagonalization of the problem from a mathematical point of view 
would be the transformation of the pCCFM equation (here using the weaker 
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constraint Q(Q — zq±)) into Mellin space: 

( Ap ccfm Af)(oo,b 2 ± ,Q 2 ) = /'^-'(A? 00 ™® Af)(x,b 2 Q 2 ) 

Jo 

JQl q ± 



x ( I q dzz"- L J ((l-z)\b ± \\q ± \)AP^ tNl (z,qi,Q 2 



+ 



I 

Jc 



*> dq]_ 

~~2 



Q Q± 



X 



r 



(33) 



Unfortunately, it turns out that the back transformation after evolution, espe- 
cially as regards the small-x region, is numerically quite problematic. Therefore, 
we take a different approach here, where both the q± and the x dependence is 
expanded into properly chosen Chebyshev polynomials. These are defined in the 
following way: 



T n (x) = cos(narccos(x)) 



(34) 



To Or) 
Ux) 
T n+1 (x) 



1 

x 



2xT n (x)-T n _ 1 (x), n>l 



r 1 T- 



dx = < 



7T 



i = j = 



(35) 



The generic advantage of the Chebyshev polynomials is now that there exists also 
a discrete orthogonality relation. Let 



x 



(JV) 



cos 



7T 



k 



N 



(36) 
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be the fc-th zero of T N (x) then one has as a discrete orthogonality relation < 
iV): 



N 



k=l 



' 






« E 

2 




(37) 


N 

V 


i = j = 





The central point is that for any arbitrary function f(x) in the interval [-1,1] the 
iV coefficients Cj given by: 



N 



C-i = — 

3 N 



E/(4 )^-i(4 j ) 
fc=i 

N 



= -E/ 



fc- 



i \ i 



COS 7T- 



cos 7r(j — 1) 



k g 
N 



(38) 



yield an approximation formula: 



/(*) 



A? 



E c fc T fc _i(x) 
.fc=i 



2 



(39) 



which is exact on all xjjf^, fc = 1, • • • , N. Therefore, one can always substitute the 
continuous integral expression for the isolation of the coefficients Cj by a discrete 
one: 

2 r 1 dx 



IT J-l y/l - 



fTk-i(x)f(x) 



N 



N 



E/ 

k=i 



k 



COS 71- 



N 



COS 7r(j — 1) 



(40) 



So, choosing a cutoff Qmax large enough, one can define a set of variables: 



t' 
t 



Hvl/Q 2 



0) 



HQUQl) 
o HQVQl) 1 



(41) 



Correspondingly for an x min small enough one can define a second set of variables: 

\nx 



y 
v 

T n {y') 



1 - 2 
1 - 2 



In x m i r 
In z 



In X m ; n 



T„ 1 - 2- 



ln z 



In i„ 



TO • 



(42) 
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In this way one can expand: 

Af(x,b ± ,q 2 ± ) = Afij(b±)c i T i _ 1 (t')c j T*_ 1 (x) , 



(43) 



ij 



where C\ = 1/2 and c k = 1 in all other cases. In this way the CCFM equation 
transforms to a simple set of linear equations: 



Af(x,b 2 ± ,Q 2 ) = 

=>A/y(&i) = 

aiji'j'ip 2 ^) = 



CCFM 



A/V,&i) + (AP 

f'ij 

ln(x min ) \n(Q 2 ma JQ 2 



Af)(x,b 2 ± ,Q 2 ) 



■Ci'Cji 



n- J-1 v'l T 1 

x I dy' I " : dt' 



Ti-iftdt f 1 Tj-Mdy 



r 1 TUXQdt r 1 
J-i VT^¥ J-i 



1 rt-2(l-y') — ln / max n 



-i 



x J ((l - z)\b ± \\qj_\)AP (z, t', m,.^ 1 ) 



/ In x m \ n . . 

ex p — o — \v -y) 



ln(x min ) \n(Q 2 max /Q 
N 2 



J ( N N 



x / dy 



k k=l k'=l 
1 rt k -2(l-yi) . 'Y m 7 ^ , 

1 



° j ^J ((l-z)|b ± ||q ± |) 

-i 

xAP CCFM (z,ql 4)^(0^(1 -( y '-y fe ,))j • (44) 



Here we used: 



tk = Vk = COS 7T- 



A: 



AT 



(45) 



The generalization for the gluon and quark singlet matrix valued expression is 
straightforward. The form of the master equation (44) has now the advantage 
that one has to handle only a simple two dimensional integration. For some 
contributions of the Altarelli-Parisi like type the expressions are diagonal in x 
and q\. In those cases more simple expressions are possible: 



Af(x,Q 2 ,bi) 



J2 c fe c //fcK & i) T fc-i(^)^„i(t Q 2) 
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(46) 
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fki{b 2 ± ) = . . . + ^2bkk'dwfk'i'{b 2 ± ) 



l'k 

2 



bkk' = J2 c k'jj 2 + qM 1 " x j) T k'-i(tj)T k -i(tj) , Xj = x 



(l-t,)/2 



a r- 2 t n \ f Q2j d< & rr u \ n2 n 

dw = 2^c v —Ti^{tj) J — 2 ^Ti'-ii^l) , Qj=Q 



IqI q 2 ± 2tt 
where tj is again the j-th zero of T N (t). 
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